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I. INTRODUCTION 

Impurity effects, arising from naturally occurring or in- 
tentionally introduced defects, have been exploited often 
as probes of the low temperature physics of various con- 
densed matter systems. The characteristic response pat- 
tern that results from the presence of impurities usually 
contains spatially resolved information which provides 
clues about the nature of the low temperature phase, 
the presence of long-range order, and the role of quan- 
tum fluctuations and other competing orders that are 
suppressed in the bulk but may be realized close to the 
impurity. Such clues are crucial in the context of mate- 
rials with strong electronic correlations whose physics is 
not as well understood as ordinary metals, and such an 
approach has therefore been heavily used in experiments 
on unconventional superconductors and Mott insulating 
antiferromagnets^ . 

On the theoretical front, the very nature of these 
strongly correlated systems makes detailed analytical so- 
lutions difficult even for simplified models. Furthermore, 
mean-field methods often fail to capture the strong cor- 
relations well enough to be reliable. In cases where ef- 
ficient numerical computations are possible, large-scale 
high-precision numerical studies can be used to access 
the physics of such systems. Other approaches that have 
had some degree of success include 'large-N' approxi- 
mation schemes, in which the symmetry group of the 
physical problem (for instance, SU(2) of spin-rotation 
invariance) is enlarged (for instance replacing SU(2) by 
SU(N)), and the resulting generalization solved by ap- 
proximations that become exact in the large- A/" limit. 
This approach has been particularly useful and influential 
in the study of quantum magnetism, and the generaliza- 
tion of the SU(2) spin rotational symmetry to SU(N) has 
played a key role in our understanding of the magnetism 
of strongly correlated Mott insulators. 

For instance, the analysis of the SU(N) Heisenberg 



model by Read and Sachde\^ ^ provided the first insights 
regarding the manner in which a Neel-ordered SU(2) an- 
tiferromagnet can be destabilized by quantum fluctua- 
tions in the antiferromagnetic order to form a Valence 
Bond Solid (VBS). In their analysis. Read and Sachdev 
considered a SU(N) generalization in which A-sublattice 
sites of the square lattice carry an SU(N) representation 
with m = 1 rows and ric columns, while 5-sublattice 
sites carried a representation with N — m = N — 1 rows 
and Tie columns. Assuming short-ranged Neel order, but 
making no assumptions about the presence of long-range 
Neel order. Read and Sachdev argued that the physics 
of such magnets can be described by a CP^~^ field the- 
ory coupled to a (compact) U(l) gauge field in two spa- 
tial dimensions, with a dimensionless coupling g that 
is a function of both Uc and the microscopic details of 
the interactions between spins (for instance, the strength 
of additional exchange or ring-exchange terms). This 
QpAT-i ^i^QQYj admits point-like space-time 'hedgehog' 
defects which correspond to monopoles in the compact 
U(l) gauge field. These hedgehogs carry Berry phases 
that depend on their spatial locations and the value of 
Tic. In the antiferromagnetic phase stabilized at small 
these hedgehogs are energetically suppressed, and their 
Berry phases can be ignored. Conversely, the VBS phase 
corresponds to a 'condensate' of these hedgehogs, with 
their position-dependent Berry phases being responsible 
for the spontaneous breaking of spatial translational sym- 
metry characteristic of the VBS phase. 

In an insightful papei^ Senthil et al. used this picture 
of the two phases to argue that the quantum phase tran- 
sition between them in two dimensional SU(2) antiferro- 
magnets can be generically a continuous transition. This 
is in conflict with expectations from the usual Landau 
theory approach, which predicts that a direct transition 
between these two phases should generically be first or- 
der. In the analysis of Senthil et al, this transition can be 
generically continuous if quadrupled hedgehog defect^, 
which are known to be irrelevant in a renormalization 
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group sense in the large- A/" limit, remain irrelevant in the 
N = 2 case. As a result of the irrelevence of hedgehog de- 
fects, the 'natural' description of this transition is then in 
terms of 'deconfined' spinon fields of a non-compact CP^ 
model (NCCP^) coupled to a non-compact U{1) gauge 
field, rather than the antiferromagnetic vector order pa- 
rameter field corresponding to the Neel ordering on one 
side of the transition or the complex order parameter 
corresponding to the VBS order on the other. 

This prediction of a 'deconfined critical point' separat- 
ing the Neel and VBS phases motivated several numerical 
studies that tried to test its validity in two quite different 
ways. One of these involved the numerical study of spe- 
cific lattice discretizations of the NCCP^ field theory it- 
self. Using such an approach, Kuklov et al. have argued 
that certain lattice discretizations of the NCCP^ field 
theory generically admit a weakly first-order transition^ 
with finite correlation length, rather than a critical point 
with diverging correlation length. However, Motrunich 
and Vishwanath have demonstrated^ that other lattice 
discretizations of the NCCP^ field theory do admit a 
critical point whose properties are expected to provide a 
good description of the generically continuous Neel- VBS 
transition. 

The second set of studies focused instead on micro- 
scopic spin models capable of displaying such Neel- VBS 
transitions. The most studied of these are the ' J — Q 
models' constructed by Sandvik^, in which a direct tran- 
sition from Neel to VBS order is driven by the competi- 
tion between four-spin plaquette interactions (Q terms) 
that favour VBS order on the square lattice and two- 
spin exchange interactions (J terms) that favour two- 
sublattice Neel order on the square lattice. While ini- 
tial studies of the J — Q model^^^ did not fully resolve 
the nature of this transition, more recent works^^^^^ have 
convincingly demonstrated the generically continuous na- 
ture of the Neel- VBS transition in this class of SU(2) 
symmetric S = 1/2 models on the square lattice. 

The most recent studies of the bulk Neel- VBS 
transitioiP^ as well as of impurity effects at this 
transitiorP^ both lead to another remarkable conclusion: 
although the Neel- VBS transition in SU(2) symmetric 
S = 1/2 systems on the square lattice is indeed con- 
tinuous, logarithmic violations of scaling behaviour, of 
the type that can arise when the critical theory contains 
marginally irrelevant operators, are also present at the 
critical point. These logarithmic effects could be the 
cause of the confiicting conclusions from the earlier nu- 
merical studies of this transition in the J — Q models, 
with some studies^ ^ interpreting their results in terms of 
a continuous transition with conventional scaling proper- 
ties, and other studies interpreting very similar data in 
terms of a weakly first-order transition^^. 

It is therefore natural to ask if this refiects the fact 
that quadrupled monopole operators, which are known 
to be strongly irrelevant in the large- limit^, actually 
become marginally irrelevant in the physical N = 2 the- 
ory. While this question cannot be directly addressed 



by computational studies in the absence of fresh insights 
on the field theory side, it does provide a strong motiva- 
tion to study similar Neel- VBS transitions on the square 
lattice in microscopic SU(N) spin models with enlarged 
symmetry {N > 2). If these transitions are continuous, 
it is of great interest to see whether they obey predic- 
tions from standard scaling theory without logarithmic 
violations. This is the question we address in this work. 

In order to do this, we focus on the SU(3) J — Q 
modePI' which has antiferromagnetic exchange interac- 
tions J between nearest neighbours and a four-spin cou- 
pling Q between sites belonging to a plaquette. In 
the language of Read and Sachde\0, the spins in this 
model carry a representation with Uc = 1 columns and 
m = 1 rows on the A-sublattice and Uc = 1 columns and 
N — m = N — 1 rows on the B-sublattice, with = 3 in 
the case of SU(3). This model generalizes the SU(N) an- 
tiferromagnetic model with purely nearest neighbour ex- 
change that was studied in earlier works^^^, who found 
that the ground state is Neel ordered for = 3 and 
A^ = 4. From the work of Lou et. al}^ and Kaul^^, it is 
known that this SU(3) model has a continuous transition 
between this Neel ordered state and a VBS ordered sin- 
glet ground state as Q is increased beyond Qc/J^ 0.502. 

Here, we study the zero temperature impurity spin 
texture that is induced by a single vacancy in the lat- 
tice at this SU(3) transition point, and confront our 
computed results with predictions from finite-size scal- 
ing theor};^^ that rely on the scale-invariance of the 
underlying quantum critical point. By a careful analysis 
of the scaling properties of the Fourier components of the 
spin texture near wavevectors k = 0, and k = Q (where 
Q = (tt, tt) is the antiferromagnetic ordering wavevec- 
tor on the square lattice), we conclude that the Neel- 
VBS transition in such SU(3) magnets on the square lat- 
tice obeys finite-size scaling forms expected at a scale- 
invariant quantum critical point, albeit with noticeable 
subleading finite-size corrections to scaling. In other 
words, the hitherto unidentified marginal operator in the 
N = 2 critical theory is no longer marginal in the A^ = 3 
case. The recent work of Kaul^^ adresses similar issues, 
albeit considering bulk properties of the transition in the 
SU(3) and SU(4) J-Q models. 

The plan of the paper is as follows. We briefiy intro- 
duce the SU(N) J — Q model in Sec.jllj We then describe 
the projector quantum Monte Carlo method used to sim- 
ulate the SU(N) J — Q model with a single non-magnetic 
impurity in Sec. |III[ first describing our choice of basis 
in Sec. |III A[ and then outlining how existing algorithms 
must be modified to account for both the existence of 
the SU(N) symmetry and of a spinless impurity. We 
then present and analyze our numerical results for the 
impurity-induced spin texture at the critical point of the 
SU(3) J — Q model in Sec. IV Finally, we discuss im- 



plications of our results and present our conclusions in 
Sec.[Vl 
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II. MODEL 

We consider a system made of Na A-sublattice sites 
with 'spins' that carry the defining representation N of 
SU(N) and Nb 5-sublattice sites at which 'spins' carry 
the conjugate representation N*. Fohowing Refs. \TT \ \Tb \ 
and [161 we write the Hamiltonian as: 

H = -J^P^J-Q P^jPku (1) 

{ij) {ij){kl) 

where Pij = —j^A^{i)A^{j), A^^ the generators of 
the SU(N) algebra satisfy the relations [A^{i)^ A^{j)] = 

^ij (^^auA^{i) — Sjj,i3A^{i)^ . The J term acts on links (ij) 

connecting nearest neighbour pairs of sites of the square 
lattice while the Q term acts on a pair of parallel near- 
est neighbour links {ij){kl) on the same plaquette of the 
square lattice. As indicated above, we choose to represent 
these generators in two different representations on the 
two sublattices of the square lattice. On the A-sublattice, 
we use the N x N matrices of the fundamental represen- 
tation N acting on states \a) {a = 1, 2, ... A^) that live 
on each yl-sublattice site, while on the 5-sublattice, the 
SU(N) generators are represented in terms of matrices of 
the complex conjugate N*of the fundamental representa- 
tion that act on the basis states \a) on each 5-sublattice 
site. 

The physics of this antiferromagnetic SU(N) model 
becomes clear upon noting that when i and j belong 
to different sublattices, Pij is actually a projector onto 
the SU(N) singlet 1 contained in the decomposition of 
N(g)N*. By analogy to the SU(2) case, the J term is thus 
expected to favour long-range antiferromagnetic order in 
the SU(N) spins, while the Q term favours a 'valence- 
bond' solid in which a particular pattern of pairs of spins 
lock into SU(N) singlets, breaking the symmetry of lat- 
tice translations. 

As is well-known, this antiferromagnetic SU(N) model 
can also be thought of as in terms of SU(2) spin S = 
{N — l)/2 variables interacting via a Hamiltonian with 
enhanced SU{2S + 1) symmetry. More precisely, we 
make the correspondence \a) = \S^ = a — S — 1) for 
a = 1,2,... A/" on the A-sublattice, and { — l)^~^\a) = 
= 5' + 1 - a) for a = 1,2,... TV on the B- 
sublattice. In this language, the SU(N) singlet state 
between a A-sublattice spin and a 5-sublattice spin 
1^*-^'^ ^ Vn ^a=i 1^)* ^ t>e rewritten as = 

^mTT^'rn=-si-^r~'\S! = m)^\S^ = -m) and thus, 
the matrix elements of Pij are explicitly given by 

^_'^^mi-\-ml-2S 

{m[,m'j\Pij\mi,mj) 2S ^'^--'^j ^rn[-m'., 

(2) 



III. METHOD 
A. Choice of basis 

Our goal is to compute expectation values of local op- 
erators in the ground-state of Hamiltonian Eq. [l] for a 
system with a given (odd or even) number of sites. For 
the SU(2) invariant, nearest neighbour 5 = 1/2 antifer- 
romagnet at Q = on a finite bipartite lattice, we know 
from the Lieb-Mattis theorem^^ that the ground-state is 
a 5'tot = singlet when the total number of sites is even, 
and is doubly degenerate corresponding to 6'tot = 1/2 
when the number of spins is odd. Unfortunately, it does 
not appear possible to come up with a simple extension 
of the Lieb-Mattis argument that would lead to such a 
result for Hamiltonian Eq. [l] when N > 2, even in the 
case with Q = 0. The main obstruction is that the 
Perron-Frobenius theorem (at the heart of the original 
Lieb-Mattis proof) cannot be used for the Hamiltonian 
in the total basis since this Hamiltonian is known to 
be reducible in each sectoi^. 

Even though we have been unable to prove the cor- 
responding statement, we nevertheless expect that the 
ground-state of Eq. [l] belongs to the SU(3) singlet sector 
when the number of sites is even. When the number of A 
sublattice sites Na equals Nb + 1, where Nb is the num- 
ber of B sublattice sites, we similarly expect, by analogy 
to the Lieb-Mattis result for S = 1/2 SU(2) antiferro- 
magnets, that the ground state of this antiferromagnetic 
Hamiltonian lies in the sector corresponding to the small- 
est irrep of SU(N) contained in the decomposition of a 
tensor product of Nb + 1 copies of N and Nb copies 
of N*. Since the smallest dimensionality irrep in this 
decomposition is N, this then suggests that the ground 
state multiplet of an SU(N) antiferromagnet with one ex- 
tra A-sublattice site will carry representation N. In the 
language of the equivalent 5* = 1 system for A/" = 3, this 
translates to the statement that the ground state for a 
system with Na = Nb + 1 has total spin Stot = 1. 

We have checked both these expectations for the pure 
Heisenberg model in one dimension with DMRG for L 
up to 64, and with Green Function Monte Carlo for the 
model Eq. [l] with L = 4 (for some values of J/Q). In 
all cases studied, the ground state was found to be in 
the expected sector. We further found (for the samples 
studied) that the lowest-energy levels in all spin sectors 
of the equivalent spin S = {N — l)/2 Hamiltonian follow 
(a weaker version of) the Lieb-Mattis ordering^^: Eq{S^ 
1)>E^{S). 

For systems with an even number of sites, the SU(N) 
singlet sector is spanned by SU(N) valence bonds, where 
each A-sublattice spin forms a SU(N) singlet with a B- 
sublattice spin^^, and the efficient projection al gorith m 
of Sandvik and Evertz can be readily generalize d l -*^^ * to 
project out the ground state written in this basis and 
sample expectation values of various operators. When 
the number of A sublattice sites is larger than the num- 
ber of B sublattice sites by one, and the ground state is 
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A^-fold degenerate (and carries the fundamental represen- 
tation N), we can proceed in a manner entirely analogous 
to the approach followed in Ref[2Tl There, the SU(2) sin- 
glet sector algorithm of Sandvik and Evertz was general- 
ized to cases where the ground state of a SU(2) symmetric 
S = 1/2 antiferromagnet has total spin 6'tot = 1/2 due 
to the presence of one extra A-sublattice site. 

To do this, we use a modified basis that leaves one A- 
sublattice site a/ 'free' (in an arbitrary state \a)af)^ and 
allows the remaining Na — 1 A-sublattice spins to form 
SU(N) singlets with the Nb = Na — 1 B-sublattice spins. 
The full basis is thus obtained by allowing all choices of 
free site a/, all choices of state \a)af for this free site, 
and all valence bond pairings possible between the other 
Na — 1 A-sublattice spins and the Nb 5-sublattice spins. 
Using an argument completely analogous to those given 
in Ref [211 for the SU(2) case, it is easy to see that this 
basis spans the ground-state sector of this SU(N) magnet 
with Nb = Na — 1 5-sublattice sites and Na A-sublattice 
sites because states from this basis can be used to con- 
struct all states of the SU(N) valence bond basis of Beach 
et. a/.l^and Lou et. a/.Elfor the singlet sector of the aux- 
iliary system obtained by adding to the odd-site system 
one extra 5-sublattice spin that carries representation 
N*. 



Ref 1211 we also need to keep track of the open string, 
since it necessarily crosses the central space- time slice, 
and therefore gives rise to an open-string in the valence 
bond loop representation of this overlap. 

Generalizing the computations of Ref. 21 to the SU(N) 
case, one can obtain modified estimators for various 
physical quantities in terms of this valence-bond loop 
representation of the overlap between {—H)'^\iIjs) and 
ii^sli—H)^. As we focus here on the impurity-induced 
spin texture, we are interested mainly in the computation 
of the expectation value {S^{t)) in the ground state with 
total z component of spin S^^^ = S. It is quite easy to 
see that this takes on a very simple expression in terms 
of the probability, Popen(r), that the open string passes 
through the site r: 

(S^r)) = (-l)"-5Pope„(r) 

This is because the contributions of all closed loops in 
the valence-bond overlap cancel out since a site covered 
by a closed loop is equally like to be in any state. 

Below, we use this to study the impurity induced spin- 
texture created by a missing-site defect at the critical 
point of the SU(3) J-Q model Eq. [l] 



IV. SPIN TEXTURE AT CRITICAL POINT 



B. Modified algorithm and estimators 



The SU(N) generalizatiorF^^ of the Sandvik-Evertz 
algorithm^^ can be modified to account for an odd num- 
ber of sites and a A^-fold degenerate ground state in a 
very simple way, following Ref [211 In this approach, the 
ground-state is obtained by stochastically sampling the 
action of {—H)^ (with some large power m ^ for a 
Lx L system of linear size L) on an initial trial state \tjjs) 
in the ground state sector with Stot = S^^^ = S with 
S = (A^ — l)/2 (in the equivalent spin representation}. 
More specifically, the 'space-time loop' representation^ 
of {iIjs\{—H)'^^\iIjs) has, apart from the closed s pace-t ime 
loops that can take on any one of 25 + 1 stateJ^^^, ex- 
actly one extra 'open-string' (as in Ref. 21) with alternat- 
ing A-sublattice and 5-sublattice nodes which are fixed 
to be in the = -\-S and = —S state respectively. 
Unlike the closed loops, the open string thus has only 
one allowed state. We use loop updates to flip between 
the 2S -\- 1 allowed states of each closed loop, exactly as 
in the SU(N) algorithm of Ref. [Til but do not change 
the state of the open string. All other aspects of the 
algorithm remain essentially unchanged, apart from the 
fact that valence-bond updates are designed in our case 
to allow the 'free site' a/ to move around with weights 
derived from the trial state \ips) written in our chosen ba- 
sis. To obtain a convenient representation of the overlap 
between (-i^)^!^^) and {ipsli-H)"^ , we follow Sand- 
vik and Evertz and keep track of space-time loops that 
cross the central space-time 'slice' to obtain the valence- 
bond loop representation of this overlap. However, as in 



A. Definition and measurements of the spin 
texture 



The impurity spin-texture {S^{t)) in the ground state 
with S^Q^ = S has an alternating part <l>^(r) that oscil- 
lates in sign between the two sublattices of the square lat- 
tice, and a uniform part <l>^(r) that varies slowly in space. 
It is therefore useful to consider the Fourier transform 

5^(k)=^(5^(r))exp(ik-r), 



defined for wavevectors k = 27rm/L (where m = 
{mx^rriy), with rrix/y = 0, 1,...L — 1). This two- 
component description of the real space texture implies 
that 6'^(k) is expected to have two peaks, one at k = 
of magnitude S (reflecting the fact that S^^^ = -\-S in 
the ground state under consideration), and the other at 
k = Q = (tt, tt) reflecting the tendency to Neel order. 

If the underlying quantum critical point is scale- 
invariant, one expects that both the uniform and the 
alternating parts of the real-space spin texture obey the 
predictions of standard finite-size scaling theory which 
assumes that the system size L is the only length- 
scale l eft in the problem at criticality. One therefore 
expect 
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FIG. 1. m = qL/27r dependence of |*S'^(k)| at the Neel-VBS transition in the SU(3) symmetric J — Q model, for k = Q + q 
(right panel) and k = q (left panel). The value of the bulk exponent 77 = 0.18 used here is obtained from a fit of the L 
dependence of |5'^(Q)p^. Insets zoom in on the small deviations from perfect scaling. 



where rj is the bulk anomalous exponent associated with 
the Neel order parameter, and are the scaling forms 
for the uniform and alternating signals, and this scaling 
form is expected to work at all |r| 1. Since it is difficult 
to unambiguously identify the alternating and uniform 
parts of the texture in real space, we find it convenient 
to go to Fourier space and note that these scaling ex- 
pectations imply that (k) obeys corresponding scaling 
laws in the vicinity of these two peaks: 

S'iq)=goim) for |q|< 7r/2 
S^Q + q) = i(3-^)/2^^(^) |q| ^ , (3) 

where m = qL/27r. 

Our strategy is thus straightforward. We use the pro- 
jector quantum Monte Carlo technique to obtain high- 
precision data on the impurity spin texture induced by a 
single vacancy in L x L samples with L even and periodic 
boundary conditions. We use the latest estimate^^ of the 
transition point J/Q\^ ~ 1.9908 (we have checked that 
our analysis holds for other slightly different estimates of 
J/Q\^). We analyze the corresponding Fourier transform 
5'^(k) near the peaks at k = and k = Q and ask if the 
scaling expectations outlined above are met. 

We begin our analysis with the peak at zero wavevec- 
tor, since this part of the analysis does not even need 
knowledge of the bulk anomalous exponent 77. We focus 
on |5'^(q)| and consider its dependence on m = qL/27r 
for various L ranging from L = 48toL = 128. The hy- 
pothesis of scale-invariance requires that all these data 
should define a single universal function of m. In other 
words, at fixed m, data from various system sizes L 
should fall on top of each other. Our results are shown 
in the left panel of Fig [l] which plots |5'^(q)| as a func- 
tion of |m| (the data is averaged over all available m that 
correspond to a given |m|). 

We see that the numerical data does seem to scale rea- 
sonably well. However, when we zoom in (inset to Fig[T]), 
small deviations from perfect scaling are visible. Next 
we investigate the behaviour near wavevector Q, using 
the best-fit value 7/ = 0.18 obtained by fitting |S'^(Q)|^^, 



and ask if |5'^(Q + 27rm/L)| collapses onto a single curve 
when all the data for various sizes L is plotted against 
m. This is shown in the right panel of Fig [l] Again, we 
see that the data does scale reasonably well, but devi- 
ations from perfect scaling are visible upon zooming in 
(inset to Fig.[T]). For both these peaks, we also note that 
the deviations from scaling are noticeably smaller than 
in the corresponding peaks for the SU(2) case (see Fig. 
1 of Ref . IT4|) . However, a detailed analysis of these small 
deviations from scaling is needed to explore this in more 
quantitative terms. 

B. Finite-size effects in the spin texture 

It is important to understand these deviations from 
scaling, and ask if they can be understood in terms of 
subleading finite-size corrections to scaling of the stan- 
dard form expected at a scale invariant quantum critical 
point. In order to do this, we note that scaling the- 
ory would predict that the dimensionless argument m 
should acquire sub-leading corrections that transform it 
to (1 + So{m){lo/ L)P^)m^ where the power po > is a 
property of the renormalization group flows towards the 
fixed point describing this scale-invariant critical point, 
So{m) is a function of the dimensionless variable m, and 
Iq is a microscopic length scale that depends on the de- 
tails of the Hamiltonian under study. A similar ansatz 
would be expected to hold for the peak centered at the 
antiferromagnetic wavevector, with ^0 replaced by ^q, Iq 
replaced by Iq and po replaced by pq. In other words, 
we modify Eq. [3] to now read: 

5-(q)=5o(m(l + C<^o(m)/LP°)) 
S^iQ + q) = L"-^gQ{m{l + F^5Q{m)/LP<^)) , (4) 

for |q| <C 7r/2. Any such ansatz that incorporates sub- 
leading finite-size effects must satisfy three conditions: 
First, the data for various sizes, when plotted against 
the scaling argument suggested by such an ansatz, must 
appear to collapse onto a single curve that defines the 




FIG. 2. Left panel: Scaling collapse of |5'^(k)| as a function of m = qL/27r at the Neel-VBS transition in the SU(3) symmetric 
J — Q model for various L in accordance with the ansatz Eqn[4]for k = Q + q (right panels) and k = q (left panels). Also 
shown are points representing Am, which fall right on top of the scaling functions qq and defined by the scaling collapse. 
We also display the |m| dependence of Bm and compare it against the expected dependence required by self-consistency. Data 
in the vicinity of both peaks are shown in linear scale (bottom panels) as well as logarithmic scale (top panels) . Straight lines 
on the log-log plots display the power-law nature of the scaling functions, with both scaling functions controlled by the single 
impurity scaling exponent rf' as indicated in the text. The linear scale of the bottom panels highlights the quality of collapse. 



corresponding scaling function g. Second, the L depen- 
dence of at each value of m should fit well to the 
form Am + B^/LP in which p is fixed by the choice made 
for the scaling variable earlier. Third, Am and 5^ must 
be consistent with expectations derived from the choice 
of scaling variable. More explicitly, A^ = ^(m), and 
Bm = I6{m)m • \/m9' These three conditions hold at 
both peaks, and the subscripts indicating the peak in 
question can be temporarily omitted from ^, J, / and p. 

Before we use this framework to analyze our data, two 
points need to be made. The first is that the range of 
sizes available is not large enough to allow a unique de- 
termination of the best fit values for A^, and p from 
fits to the L dependence at each m. This is because the 
deviations from perfect scaling are rather small, and as a 
result, one needs a very large range in L to pin down the 
best fit p. With our limited range in L, a range of p gives 
equivalently good fits for the L dependence at each m. 
Each choice of p in this range fixes a best-fit value of A^ 
and Bm when one fits the L dependence at individual m. 

Likewise, the deviations from perfect scaling are too 
small to uniquely fix p, S and / by demanding the best 
scaling collapse for all the data plotted against the scaling 
variable corresponding to these choices for p, 5 and One 



obtains scaling collapse that looks equally good to the eye 
for a range of p, and for at least two simple choices for 
the function 6: ^(m) = 1 and (5(m) = |m|. The best one 
can do is to fix / uniquely after choosing a p within this 
range and making one of these two choices for 5. 

In view of this uncertainty, the self-consistency condi- 
tions outlined above (which relate the best-fit Am and 
5rn to the choices made for 5 and / which define the cor- 
responding scaling variable, and to the function g defined 
by the resulting scaling collapse) are particularly useful 
in fixing the form of the finite-size corrections from the 
available data. We have therefore systematically scanned 
over the range of p that yields a viable description for 
both simple choices of the function S. For each choice of 
S and p, we have determined the value of / that yields the 
best scaling collapse, as well as obtained the best-fit Am 
and 5m from fits to the L dependence at each individual 
m. We have then tested for self-consistency as outlined 
above. 

Based on this extensive study, we conclude that fi- 
nite size corrections to perfect scaling take on the form 
corresponding to So{m) = SQ{m) = |m|, Iq = 3.8(1), 
Iq = 6.5(1), po = 1.0(2), and pq = 1.2(1). Furthermore, 
the scaling functions ^o(^) and ^q(x) both take on power- 
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law forms for not too small x, consistent with the expec- 
tations from the scaling theory for impurity respons^^ 
which fixes the power-law forms of both functions in 
terms of a single impurity exponent r]': go{m) ~ |m|~^ 
and gQ{m) - |m|-(3-^+^')/2. 

This scaling collapse is shown in Fig. [2j while the cor- 
responding fits to the data at fixed m, from which the 
values of and are obtained, are shown in Fig. |3] 
The Ayyi thus obtained are marked in Fig. [2] and seen 
to fall right on top of the composite curves produced by 
the data collapse which define the corresponding scal- 
ing function go and ^q. These scaling functions are fit 
to power law forms expected from the scaling theory for 
the impurity spin texture, and both scaling functions are 
seen to be well-described by power-laws determined by a 
single impurity exponent rj' = 1.0(1), as expected from 
the scaling analysis. Also shown in Fig. |2]is the compar- 
ison between the best fit 5^ obtained from the fits in 
Fig. [3] and their expected |m| dependence obtained from 
the self-consistency. Again, we see that the |m| depen- 
dence of Bm is consistent with that expected from the 
consistency condition. 

Thus, conventional finite-size corrections to scaling ap- 
pear to give a completely consistent account of our data 
for both peaks in the Fourier transform of the impurity 
spin texture. In this context, it is interesting to note that 
multiplicative corrections (which amount to violations of 
scaling) appear to be favoured over such conventional 
additive corrections to scaling by the author of Ref. [T3l 
when describing the behaviour of the spin stiffness at the 
transition in the pure system. 



of our recent stud}'^^ of the spin texture in the SU(2) 
symmetric S = 1/2 J — Q models on the square lat- 
tice. There, we found that any attempt to fit the de- 
viations from scaling to this conventional form resulted 
in unphysically large / and very small p in the SU(2) 
case. As a result, one was forced in the SU(2) case to 
consider either an unconventional form of scaling argu- 
ment proportional to m/L^ with very small p, or admit 
the possibility of logarithmic violations of scaling, as re- 
flected in an argument of the form m/log(L//) for the 
scaling function. These two possibilities are in fact nu- 
merically indistinguishable from each other in the SU(2) 
case over the range of available sizes. Since the former 
has no known basis while the latter is a known to be a 
possible consequence of marginally irrelevant operators 
in the renormalization group description of phase tran- 
sitions, the most natural conclusion consistent with the 
numerical evidence was that there are logarithmic viola- 
tions of scaling d ue to t he presence of a marginal operator 
in the SU(2) cas4 ^^^ . That this is not the case at the 
SU(3) critical point is the central message of our work. 

Our results thus support the scenario that some op- 
erator that is irrelevant at fixed points describing Neel- 
VBS transitions for N > 3 becomes marginal at or near 
N = 2. We hope these results will motivate further study 
of the field-theoretical description of the critical theories 
underlying these interesting, but still only partially un- 
derstood, quantum phase transitions. 
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